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ABSTRACT 



X 

Nucleon structure functions are studied within the chiral soliton approach to the boson- 
ized Nambu-Jona-Lasinio model. The valence quark approximation is employed which is 
justified for moderate constituent quark masses (~ 400MeV) as the contribution of the va- 
lence quark level dominates the predictions of nucleon properties. As examples the unpolar- 
ized structure functions for the up and up scattering and the structure functions entering the 
Gottfried sum rule are discussed. For the latter the model prediction is found to reasonably 
well agree with a corresponding low-scale parametrization of the empirical data. 
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1. Introduction 

It has been a long standing puzzle how to establish the connection between the chiral 
soliton picture of the baryon, which essentially views baryons as mesonic lumps, and the 
quark parton model, which regards baryons as composites of almost non-interacting, point- 
like quarks. While the former has been quite successful in describing static properties of the 
nucleon, the latter, being firmly established within the context of deep inelastic scattering 
(DIS), has been quite successful in predicting the spin average DIS nucleon structure func- 
tions. The apparent difference between models for the nucleon like the bag model QXJ] , which 
have previously been employed to study structure functions and soliton models is 

the fact that in the latter the nucleon wave-function only appears as a collectively excited 
(topologically) non-trivial meson configuration rather than as a product of Dirac spinors. 
In this letter we calculate structure functions in the Nambu-Jona-Lasinio (NJL) [H| chiral 



soliton model [pi] , |12jj where the hadronic currents are formally described in terms of quark 
degrees of freedom which themselves are functionals of the solitonic meson fields. Since the 
present study is the first step towards computing nucleon structure functions from a chiral 
soliton we will adopt a simplifying valence quark type of approximation (to be defined after 
eq (^)) and leave a more complete exploration to future studies. 

As in the original study of structure functions for localized field configurations, the 
structure functions are most easily accessible when the current operator is at most quadratic 
in the fundamental fields and the propagation of the interpolating field can be regarded as 
free. Although the latter approximation is well justified in the Bjorken limit the former 
condition is difficult to satisfy in soliton models built from mesonic fields. In such models 
the soliton is a non-perturbative object involving all orders of the fundamental pion field. 
Hence the current operator is not confined to quadratic order. In models where mesons are 
fundamental fields {e.g. the Skyrme model |L3], [TJ]], the chiral quark model of ref. [15| or the 



chiral bag model flqln) structure functions are exceedingly difficult to obtain. In this respect 
the chirally invariant NJL model is advantageous because it is entirely defined in terms of 
quark degrees of freedom. This makes the evaluation of the required commutator (see eq 
(|3]) below) feasible. Nevertheless the quark currents become uniquely (up to regularization) 
defined functionals of the meson fields. The Lagrangian of the NJL model reads 

£ = q(i$ - m°)q + 2G NJL £ ([q^qf + (q^q) 2 ^ . (1) 

Here q, rh° and Gnjl denote the quark field, the current quark mass and a dimensionful 
coupling constant, respectively. Functional bosonization jT7| yields the action 



1 f 

A = Tr A log(iD) + — / d 4 x tr (m° (M + M f ) - MM" 1 ") , (2) 

D = i$ - (M + M f ) - 75 (M + M f ) . (3) 



The composite scalar (S) and pseudoscalar (P) meson fields are contained in M = S + iP, 
and appear as quark-antiquark bound states. For regularization, which is indicated by the 

*In the cloudy bag model the contribution of the pions to structure functions has been treated perturba- 
tively|J. 
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cut-off A, we will adopt the proper-time scheme |18|]. The free parameters of the model are 
the current quark mass m°, the coupling constant Gnjl and the cut-off A. When expanding 
A to quadratic order in M these parameters are related to the pion mass, m n = 135MeV and 
decay constant, f n = 93MeV. This leaves one undetermined parameter which we choose to be 
the vacuum expectation value m = (M). For apparent reasons m is called the constituent 
quark mass. It is related to m°, Gnjl and A via the gap-equation, i.e. the equation of 
motion for the scalar field S. The occurrence of this vacuum expectation value reflects the 
spontaneous breaking of chiral symmetry. 

We will approach the computation of structure functions in the NJL model by first 
briefly reviewing the kinematics of the Bjorken limit and the NJL soliton in sections 2 
and 3, respectively. In section 4 we will work out the valence quark approximation to the 
unpolarized structure functions. The numerical results will be presented in section 5. Finally 
section 6 not only serves to summarize these studies but also to propose further explorations. 



2. Kinematics 

The starting point for computing nucleon structure functions is the hadronic tensor 



W ab (q) = — 



e^<(N\ tofl.J?(0)l \N). 



(4) 



Here \N) refers to the nucleon state and J%(0 = <?(07/Ai<?(0 is the hadronic vector current. 
In the context of weak interactions we take J^(£) = g(£)7 M (l — J5)t a q(£). We denote by 
t a (a = 0, ..,3) the flavor operators, which have to be chosen appropriately for the process 
under consider at ionP]. The unpolarized structure functions are related to the symmetric 
piece, Wffl = (W^ u + W vft )/2, which is parametrized by two scalar form factors, 



W<?(?) 



-9^ + ^rl Wi(xBj,Q 2 ) + (P, - q, 1 -^] [P v - q v ^-f] W 2 (x Bi , Q 2 ) . (5) 



Here P M refers to the nucleon momentum and Q 2 = —q 2 . Furthermore quantities suitable 
to study the Bjorken scaling have been introduced: v = P ■ q/M]y and xb] = Q 2 /2M^v. 
Introducing the projection operators pf 



Af = - 
1 2 



-g 



V 



M 2 N 



Ar 



-g 



hi/ 



3 7) 



pV- pv 



(6) 



with i] = 2M N uxBj/{2M N i'XBj + v) enables one to straightforwardly extract the form factors 



W t (x B] ,Q 2 )=ArW^(q) , i=l,2 



(7) 



When discussing Bjorken scaling a slightly different definition of the form factors 



FxixBj, Q 2 ) = M N Wi(xBj,Q 2 ) and F 2 (x Bj , Q 2 ) = vW 2 (x Bi , Q 2 



(8) 



is commonly considered. The Bjorken limit corresponds to the kinematical regime 



\q\ — MnXb] with \q\ — > oo . 



(9) 



^For simplicity we omit the flavor index when not relevant. 
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Finally the structure functions are obtained as the Bjorken limit of the form factors (Q) 

Fi(xBj) = Hm Fi{x B] ,Q 2 ) , i = 1,2 . (10) 

3. The Nucleon State in the NJL Model 

As the NJL model soliton has exhaustively been discussed in a recent review article 



lzf we only present those features, which are relevant for the computation of the structure 
functions. 

The chiral soliton is given by the hedgehog configuration of the meson fields 

Mn(x) = m exp (it ■ xQ(r)) . (11) 

In order to compute the functional trace in eq (|2]) for this static configuration a Hamilton 
operator, h is extracted from the Dirac operator (^). That is, D = i n fo(dt — h) with 

h = a ■ p + m exp (275 r • xQ(r)) . (12) 

We denote the eigenvalues and eigenf unctions of h by e M and respectively. In the proper 



time regularization scheme the NJL model energy functional is found to be [JTTJ, [L2| 

E[0] = ^ (1 + sgn(e v )) + ^ T 2 ~^= 3 E ™P ("^) + <fl I ^ (1 - ©(0) - (13) 

with iVc = 3 being the number of color degrees of freedom. The subscript "v" denotes the 
valence quark level. This state is the distinct level bound in the soliton background with 
— m < e v < m. The chiral angle, G(r), is obtained by self-consistently extremizing E[Q] 
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Nucleon states possessing good spin and isospin quantum numbers are generated from 
the soliton by taking the zero-modes to be time dependent |fL4 



M(x,t) = A(t)M H (x)A\t) , (14) 
which introduces the collective coordinates A(t) e SU(2). The action functional is expanded 



[TTJ in the angular velocities 

2A\t)A(t) = it • n. (15) 

In particular the valence quark wave-function receives a first order perturbation 

* v (aj, t) = e-^A(t) I * v (aj) + W ^(g) ^ |T ^" |v) 1 =: A(t)^(x) . (16) 

Here ip v (x) refers to the spatial part of the body-fixed valence quark wave-function with the 
rotational corrections included. Nucleon states \N) are obtained by canonically quantizing 
the collective coordinates, A(t). By construction these states live in the Hilbert space of a 
rigid rotator. The eigenf unctions are Wigner D-functions 

(A\N) = ^D)i% 3 (A) , (17) 
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with J3 and J3 being respectively the isospin and spin projection quantum numbers of the 
nucleon. 



4- Unpolarized Structure Functions in the Valence Quark Approximation 

The starting point for the computation of the unpolarized structure functions is the 
hadronic tensor in the form suitable for localized fields |2| 

-— 4 S, pva k p sgn (k ) 6 (k 2 ) J ^ dt e'^ 4 ** 
x J d 3 xi J d 3 x 2 exp [—i(k + q) ■ (x% — x 2 )] 

x (N\ t)tit m -f^(x 2 , 0) - V(x 2 , 0)t m t n ^{ Xl , t)} \N). (18) 

Here S^ pva = g w g va + g^g vp - 9^g P a denotes the symmetric combination of 7^7^ and 
( = 1(2) for the structure functions associated with the vector (weak) current. As explained 
in the preceding section the matrix element between the nucleon states (\N)) is to be taken 
in the space of the collective coordinates. In deriving the expression fll8|) the free correlation 
function for the intermediate quark fields has been assumed 0. This reduces the commutator 
[J^Xijt), Jl(x 2 ,0)} of the quark currents in the definition (||) to objects which are merely 
bilinear in the quark fields. In the Bjorken limit @ the momentum, k, of the intermediate 
quark state is highly off-shell and hence not sensitive to momenta typical for the soliton 
configuration. Thus the use of the free correlation function is a good approximation in this 
kinematical regime. Accordingly, the intermediate quark states are taken to be massless, cf. 
eq (|T%D . In the next step the form factors Wi are extracted according to eq (0) . Noting that 

Sppva-^gpa and Aff S upva—^rjg pa the Callan-Gross relation follows immediately, i.e. 
-^(^Bj) = 2a/Bji*i(£Bj)- It thus suffices to only consider the structure function Fx(xBj). 

Since the NJL model is formulated in terms of quark degrees of freedom, quark bilinears 
as in eq ( |18|) can be computed from the functional 
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(q(x)Vq(y)) = J DqDq q(x)Tq(y) exp J d 4 x' C 

■ J DqDq exp (i J d'x'd^y' [6\x - y)C + a(x' , y')q{x')Yq{y')] ) \ a ^ y)=Q , (19) 



iSa(x, y) 

where T is a suitable Dirac and/or isospin matrix. The introduction of the bilocal source 
a(x, y) facilitates the functional bosonization upon which eq ( [I9|) takes the form 

d Tr A log (5\x - y)D + a(x, y)T)) , . (20) 



5a(x,y) 



a(x,y)=0 



The operator D is defined in eq (S). From this discussion it is obvious that structure functions 
are most easily obtained within models which can completely be formulated in terms of quark 
degrees of freedom where the form of the current operator is not altered by the interactions^]. 

* Otherwise matrix elements of operators have to be computed, which are more complicated than the 
bilocal quark bilinear q(x)Tq(y). 
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The correlation (q(x)Tq(y)) depends on the angle between x and y. Since in general 
the functional (|1^) involves quark states of all angular momenta (I) a technical difficulty 
arises because the angular dependence has to be treated numerically. The major purpose 
of the present letter is to demonstrate that structure functions can be computed from a 
chiral soliton. With this in mind we will adopt the valence quark approximation where only 
quark orbital angular momenta up to I = 2 are relevant. From a physical point of view 
this approximation is justified at least for small constituent quark masses because in that 
parameter region the nucleon properties are dominated by the valence quark contribution 
I2" |. We define the valence quark approximation to the structure functions by restrict- 



ing the quark configurations in (^) to the iso-rotating valence quark wave-function (|l6|) , 
accordingly substituting the valence quark wave-function (|16|) into eq (|18|). 

When extracting the structure function Fi(xbj) the integrals over the time coordinate 
can readily be done yielding the conservation of energy for forward and backward moving 
intermediate quarks. Carrying out the integrals over k and k = \k\ yields for the isovector 
part of the structure function 

F( =1 (x Bi ) = (Nc^r(N\D 3i J dn k k 2 ^l(p) (l + a-k) r^ v (p)\ ^^tr (r 3 ^ m ) 

-$(-p) (l + a-k) 7#v(-p) h tr (r 3 t m ti) } \N) , (21) 

V / k=qo— e v I 

where p = k + q. Nc appears as a multiplicative factor because the functional trace (|20|) 
includes the color trace as well. Furthermore the Fourier transform of the valence quark 
wave-function 

Vv(p) = J-^~ ^vO) exp (ip ■ x) (22) 

has been introduced. The isoscalar part of the structure function, F/ =0 (a;Bj), is straightfor- 
wardly obtained from eq ( f21~l) by replacing with unity and omitting the isospin matrices 
Tj. Note that the wave-function ip v contains an implicit dependence on the collective coor- 
dinates through the angular velocity fl, cf. eq (|l6l). 

The dependence of the wave-function i/j(+p) on the integration variable k is only implicit. 
In order to carry out this integration it is most convenient to choose the external momentum 
along the z-axis, i.e. q = qz. In the Bjorken limit the integration variables may then be 
changed to [0] 

k 2 dQfc = pdp d$ , p = |p| , (23) 

where $ denotes the azimuth-angle between q and p. The lower bound for the p-integral is 
adopted when k and q are anti-parallel: p™ m = \M N x B j =p e v | for k = q ± e v , respectively. 
The wave-function ip(+p) acquires its dominant support for p < M^. Hence the integrand 
is different from zero only when q and k are anti-parallel and we may take k = —z. This is 
nothing but the light-cone description for computing the structure functions 0] . The valence 
quark state possesses positive parity yielding ip(—p) = joip(p). We now have arrived at the 
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final expression for the isoscalar and isovector parts of the unpolarized structure function in 
the valence quark approximation 



F I=0 



( (Fl(x Bi ) - F[(x Bi 
N c 



-^(N\ pdp d$^(p T )(lTa 3 )Vv(p T )|iV)tr[^ ro ] (24) 

An Jm n \x^\ Jo 



=F 

Mm r°° r 2n 

N c -^(N D 3i / pdp I di> 



x$(p Jr, (1 T a 3 ) Vv(p =)\N)ti 



T 3 



tlt-m 



(25) 



where x± = x B j ±e v /M N and cos(6p) = M N x±/p. The polar-angle, 0^ , between p ± and q 
is fixed for a given value of the Bjorken parameter, x B - y Hence the wave-function depends 
implicitly on x B j because Vv(p±) — Vv(p, Qp, $)• 

Turning to the evaluation of the nucleon matrix elements defined above we first note that 
the Fourier-transform of the wave-function is easily obtained because the angular parts are 
tensor spherical harmonics in both coordinate and momentum spaces. Hence, only the radial 
part requires numerical treatment. Performing straightforwardly the azimuthal integrations 
in eqs fl2"l|) and fl2"5D reveals that the isoscalar part, -F/ =0 , depends solely on the classical 
part of the valence quark wave-function, ty . Thus F[ =0 is identical for all nucleon states. 
On the other hand the isovector part, F( =1 , is linear in the angular velocity, O. Since 
the ^-direction is distinct the collective quantities appear as combinations of .033^3 and 
D 3i Qi. When quantizing the collective coordinates these combinations are substituted by 



the nucleon spin operator yielding [J4 



{N\D S3 n 3 \N) 



h 



3a 2 



and {N\DM\N) = ~4 . 



a- 



(26) 



Here ^3 = ±(1/2) is the nucleon isospin projection and a 2 refers to the moment of inertia 
of the soliton. For consistency we constrain it to the valence quark contribution, a 2 , cf. eq 



(|3l|). The isovector part is obviously proportional to the isospin projection but independent 
of the spin projection, as expected for unpolarized structure functions. It is convenient to 
define structure functions f±(x B j) with the nucleon matrix elements already computed via 



F°(x Bj ) = f° ± (x Bi ) and i^(x Bj ) = 2J 3 /i(x B j) . 



(27) 



5. Results 

In figure 1 we display the unpolarized structure functions f^ 1 for a constituent quark 
mass of m = 350MeV. In that case the valence quark contribution to the moment of inertia 
is about 86%. This shows that the vacuum is only moderately polarized and that the valence 
quark approximation is well justified. Here we assume the experimental value (940MeV) for 
the nucleon mass. We observe that the structure functions are well localized in the interval 
< x B j < 1. The result that the structure functions slightly exceed x B j = 1 is common to 
approaches which treat the nucleon as extended objects. In the context of bag type models 
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Isoscalar Structure Functions 



Isovector Structure Functions 





Figure 1: The unpolarized structure functions f^ 1 as functions of the Bjorken variable xbj- 

various projection techniques have been proposed || |J to remedy this problem. This, 
however, is not the central issue of this paper. 

The results displayed in figure 1 are the central issue of this calculation and it is of 
great interest to compare them with the available data. In this context we consider the 
structure functions for electron nucleon scattering. The associated isospin matrices are 
tah = t b t a = (5 + 3r 3 )/18 yielding 



Ff = I (5(/J - /°) - 6/ 3 (4 - fl) 



1 



(28) 



where the second equation results from the Callan-Gross relation. As for all effective low- 
energy models of the nucleon, the predicted results are at a scale lower than the experimental 
data. In order to carry out a sensible comparison either the model results have to be evolved 
upward or the QCD renormalization group equations have to be used to extract structure 
functions at a low-renormalization point. The latter procedure has been employed in ref plfl p] 
to make available a low-scale parametrization of the empirical data on F£ N . From figure 
2 we observe that the NJL model prediction for F 2 ep — F 2 en reproduces the gross features 
of this parametrization although the maximal value of the prediction is a bit too large. 
On the other hand the low-scale value are more enhanced at small xej- To illustrate the 
origin of the bumb at XBj ~ we have also included the low-scale parametrization with 
the a s -corrections omitted, cf. eq (7) of ref |^T|. These are actually the starting point for 
computing the low-scale parametrization. When including the a s -corrections the integral 
(|29|) is forced to remain unchanged. As the a 2 ~corrections shift the structure functions to 



§These authors also provide a low scale parametrization of quark distribution functions. However, these 
refer to perturbatively interacting partons. Distributions for the NJL-model constituent quarks could in 
principle be extracted from eqs. (p^)-(p5|). It is important to stress that these distributions may not be 
compared to those of ref |2M because the associated quarks fields are different in nature. 
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(Anti)Neutrino Scattering 



m = 400MeV 
m = 450MeV 
Glueck et al w/o a 5 [21] 
Glueck et al w a s [21] 




0.00 0.25 0.50 0.75 1.00 1.25 




Figure 2: The valence quark approximation to the unpolarized structure functions as func- 
tions of the Bjorken variable XBy Left panel: The prediction on the Gottfried sum for two 
values of the constituent quark mass m. We compare with the low-scale parametrization of 



ref pi. Right panel: F" p and F v x v for m = 450MeV. 



larger xb] the (artificial) bumb at XBj ~ emerges. As an aside we would like to mention 
that the agreement between the NJL-model predictions and the parametrized structure 
functions is better when the a s -corrections are omitted. This indicates that a fine-tuning 
of the low-scale momentum might improve the agreement even more. 

With regard to the vacuum contribution it should be emphasized that it will not simply 
add to the valence quark piece because when computing the isovector structure functions 
we have substituted al < a 2 in eq (EH). For m = 400(450)MeV we find al/a 2 = 78(72)%. 



We observe that for the combination F 2 — F™ the agreement with the parametrization 



21] improves as m increases. We conjecture that this feature survives when the vacuum 



contribution is included because the moment of inertia enters the denominator of F 2 P — F 2 n . 



Furthermore the integral 



S, 



G 



dx 



Bj 



(F? - F- 



Bj 



dx B] {F{ v - F{ n ) = 0.29 (0.27) 



(29) 



agrees reasonably well with the empirical value Sq = 0.235 ±0.026 122] for the Gottfried sum 
rule. In particular the deviation from the naive value (1/3) |23| is in the correct direction. 

For the weak scattering precesses up and up we demand t a tb = (1 ± r 3 )/2 yielding the 
linear combinations 



Fr = 2(f + -f +fl + f 1 _) and Ff = 2 (f° + - f_ - f + - f_ 



(30) 



which are also plotted in figure 2. Although our wave-functions (cf. section 3) are quite 
different from those in the bag model the shape of the structure functions is similar. In 
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particular the structure functions F" v,vv do not vanish at xbj = in both models. Despite 
that we essentially take only one quark eigenstate into account, we find a clear smearing of 
the structure functions. This shows that relativistic effects, i.e. a sizable lower component 
of the valence quark wave-function, play a significant role. These effects also cause the 
maximum of the structure to be shifted from e v /M pa 0.26 to about 0.37. As in the bag 
model calculation of ref we find that _Ff p is negative in the vicinity of XBj = 0. This 
appears to be linked to the omission of the vacuum states when computing the hadronic 
tensor (p~8|). 

Let us briefly comment on the Adler sum rule. Note that in the Bjorken limit, where the 
Callan-Gross relation is satisfied, the Alder and Bjorken sum rules are equivalent. It is an 
easy matter of exercise to verify that 

Thus the Alder sum rule is satisfied once we assign the moment of inertia to its valence 
quark contribution, a^. It is obvious that this sum rule will be recovered without this 
restriction when the contribution of the polarized vacuum is included in the evaluation of the 
functional trace ( p0|) . The Adler sum rule also serves as a test for our numerical treatment. It 
furthermore manifests the parton model interpretation because adopting unity as the upper 
boundary of the integral ([H]) saturates this sum rule already by 99% for the parameters 
used here. We should mention that the momentum sum rule is not satisfied in the valence 
quark approximation. The analytical proof of the momentum sum rule involves the classical 
equation of motion for the chiral field. As the polarized vacuum contributes to this equation 
it is obvious that including only the valence quark level in the calculation of the structure 
functions violates this sum rule. Numerically, however, this violation is small. For example, 
for m = 450MeV we are missing about 20%. This number decreases with the constituent 
mass and may be interpreted as the momentum carried by the polarized vacuum^} 

6. Summary and Outlook 

The present study is intended as the first step towards clarifying the connection between 
the chiral soliton picture of the nucleon and the quark parton description. This has also to 
be regarded as an attempt to combine the phenomenologically successful concept of chiral 
symmetry with the quark parton model. For this purpose we have presented a first calcu- 
lation of nucleon structure functions in the Bjorken limit from a chirally symmetric model. 
The mean field quark wave-functions in the background of a chiral soliton represent a non- 
trivial coupling of spin and isospin to the so-called grand spin. Baryon states possessing 
good spin and isospin are subsequently generated by cranking the soliton, see eq flIS|). As 
a consequence there are rotational corrections to the mean field predictions of the structure 
functions. These corrections contribute to the isovector part of the unpolarized structure 
functions and are mandatory to reproduce the Adler sum rule. This form of the nucleon 

"This consistency check requires one to use the soliton mass (~ 1.2GeV) rather than the experimental 
value for the nucleon mass. 
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wave-function constitutes a major difference to quark models which are not based on a non- 
trivial chiral field as e.g. the bag model where baryons are described as direct products of 
spin and isospin eigenstates. In order to establish the connection between the two pictures 
we have (as a first step) restricted ourselves to the valence quark approximation. The results 
are in reasonable agreement with bag model results (which does not include the Dirac sea 
either) and the empirical value of the Gottfried sum rule. Also the low-scale parametrization 
of the combination — F™, which enters this sum rule, is satisfactorily reproduced. 

This encourages further studies in various directions. It is obvious that for the full 
computation of the structure functions the polarization of the vacuum quark states has to 
be incorporated. Although for small constituent quark masses the soliton is dominated by 
the valence quark configuration it will be very interesting to have available direct access 
to the sea quarks. The regularization of the functional trace ( p0|) will be rather involved. 
Fortunately the Adler and momentum sum rules may be employed to perform consistency 
checks. This study will illuminate whether (and how) the vacuum contribution violates the 
identity /+(0) = /£(0), which is also observed in the bag model @. Within the parton 
model picture the sea quarks cause a violation of this relation. 

Of special interest are the polarized structure functions which are to be extracted from 
the anti-symmetric part of the hadronic tensor Wffl = (W^ v — W Ufl )/2. The smallness of 
the first moment of the associated flavor singlet structure function is known as the proton 
spin puzzle. Since almost all chiral soliton models provide a reasonable explanation of this 
puzzle the computation of the entire structure function will provide further understanding 
how the nucleon is built up from its constituents. 

One wonders whether the functional trace ( |20| ) has a suitable interpretation in chiral mod- 
els with mesons as the fundamental fields. Although it is possible to identify quark bilinear 
quantities in such soliton models via saturation of the Ward-identities, the analogues of the 
quark bilinears are always local. Hence models with fundamental meson fields contributing 
to the currents seem to be less tractable for calculating structure functions. In this respect 
an expansion of the functional trace (pOf) in derivatives of the chiral field might provide an 
effective operator suitable to compute structure functions in purely mesonic soliton models. 



However, first investigations along this line have given disappointing results |24fl . This might 
be related to the failure of the gradient expansion in the soliton sector. 

Gluonic effects are known to significantly contribute to the structure functions, they may 
even cause some of them to be singular^ at xb^ = 0. Such singularities will not appear in the 
soliton model calculation (neither do they in the bag model calculation ||). Hence a further 
study of the structure functions may provide some insight how to effectively incorporate 
gluonic degrees of freedom in NJL-type models. 
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